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ABSTRACT 



Aims. To understand stellar magnetism and to test the validity of the Babcock-Leighton flux transport mean field 
dynamo models with stellar activity observations 

Methods. 2-D mean field dynamo models at various rotation rates are computed with the STELEM code to study the 
sensitivity of the activity cycle period and butterfly diagram to parameter changes and are compared to observational 
data. The novelty is that these 2-D mean field dynamo models incorporate scaling laws deduced from 3-D hydrody- 
namical simulations for the influence of rotation rate on the amplitude and profile of the meridional circulation. These 
models make also use of observational scaling laws for the variation of differential rotation with rotation rate. 
Results. We find that Babcock-Leighton flux transport dynamo models are able to reproduce the change in topology of 
the magnetic field (i.e. toward being more toroidal with increasing rotation rate) but seem to have difficulty reproducing 
the cycle period vs activity period correlation observed in solar-like stars if a monolithic single cell meridional flow is 
assumed. It may however be possible to recover the P cyc vs P ro t relation with more complex meridional flows, if the 
profile changes in a particular assumed manner with rotation rate. 

Conclusions. The Babcock-Leighton flux transport dynamo model based on single cell meridional circulation does not 
reproduce the P cyc vs P ro t relation unless the amplitude of the meridional circulation is assumed to increase with 
rotation rate which seems to be in contradiction with recent results obtained with 3-D global simulations. 
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1. Introduction 

The Sun is a magnetic star that possesses a well defined 
cyclic activity period of 22 years. It is believed that the 
origin of this magnetism and regular activity is linked 
to a global scale dynamo operating in and at the base 
of the solar convection zone (Parker 1993). In order to 
model and understand this global scale dynamo, a use- 
ful approach has been to make use of the mean field dy- 
namo theory (jMoffatt 19781 Krause & Raedler 1980). This 
method has the advantage that it only deals with the 
large scale magnetic field, assuming some parameterization 
of the underlying small scale turbulence and magnetism. 
Computing those physical processes self-consistently would 
otherwise require costly full 3-D magnetohydrodynami- 
cal (MHD) simulations (Cattaneo 1999; Brun, Miesch & 
Toomre 2004). Among the various mean field dynamo mod- 
els, the Babcock-Leighton flux transport dynamo models 
have recently been the favored ones and have demonstrated 
some success at reproducing solar observations assuming ei- 
ther an advection or a diffusion dominated regime for the 
transport of field from the surface down to the tachocline 
(Dikpati & Charbonneau 1999; Nandy & Choudhuri 2002; 
Dikpati et al. 2004; Chatterjee, P. et al. 2004; Charbonneau 
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2005; Bonanno et al. 2005; Jouve & Brun 2007; Yeates et 
al. 2008). 

It is legitimate to ask if such models can be applied to 
other stars and lead to the same agreement with the ob- 
servations of for instance solar type stars that may possess 
different rotation rates and activity levels. 

1.1. Observational Constraints 

Observations of the magnetism of 'solar type' stars, i.e. 
stars possessing a deep convective envelope and a radiative 
interior (late F, G, K and early M spectral type) are becom- 
ing more and more available (Giampapa 2004) and may pro- 
vide additional constraints on our understanding of global 
scale stellar dynamos. One difficulty of such observational 
programs is that they require long term observations since 
stellar cycle periods are likely to be commensurate to the 
solar 11-yr sunspots cycle period (or 22-yr for a full cycle 
including two polarity reversals of the global poloidal field). 
Thanks to the data collected at Mount Wilson Observatory 
since the late 60's, such data is available (Wilson 1978; 
Baliunas et al. 1995). Among the sample of 111 stars (in- 
cluding the Sun as a star) originally observed between the 
F2 to M2 spectral type, it is found that about 50% of the 
stars possess a cyclic activity, with cycle (starspot) periods 
varying roughly between 5 to 25 yrs, i.e. between half to 
twice the sunspot cycle period. They further indicate that 
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among the inactive stars of the sample some are likely to 
be in a quiet phase (as was the Sun during the Maunder 
minimum). Overall activity cycles seems to be more fre- 
quent for less massive stars K than for F stars. More re- 
cent observational programs have been pursued that now 
even provide information on the field topology as a func- 
tion of the rotation rate, such as the one using the Espadons 
and Narval instruments and the Zeeman Doppler Imaging 
technique (Donati et al. 2006). Applying this observational 
technique over a sample of four solar analogues with rota- 
tion rate varying from one to three times solar, Petit et al. 
(2008) have shown that the field amplitude increases as a 
function of the star's rotation rate and, more importantly, 
becomes more and more dominated by its toroidal compo- 
nent (modulo possible bias in the observational technique 
used). If such a trend is confirmed, i.e. that the field topol- 
ogy is becoming more toroidal with increasing rotation rate, 
it is a very important and instructive result and puts strong 
constraints on the dynamo models. 

The systematic analysis of stellar magnetism data re- 
vealed that for solar type stars there is a good correlation 
between the cycle and rotation periods of the stars and 
that correlation is even stronger when using the Rossby 
number (Ro = P rot /r) that takes into account the convec- 
tion turnover time r at the base of the stellar convective 
envelope (Noyes et al. 1984; Soon et al. 1993; Hempelmann 
et al. 1995; Baliunas et al. 1996). As the star rotates faster, 
its cycle period is found to be shorter. Typically Noyes et 
al. (1984) found that P cyc oc P r n o£ , with n = 1.25 ± 0.5. 
Based on an extended stellar sample Saar & Brandenburg 
(1999) and Saar (2002) have argued that there is actually 
two branches when plotting the cycle period vs the rota- 
tion period of the stars. They make the distinction between 
the primary (starspot) cycle and Gleissberg or grand min- 
ima type modulation of the stellar activity. For the active 
branch they found an exponent n ~ 0.8 and for the inac- 
tive stars n ~ 1.15 (Charbonneau & Saar 2001). It is also 
found that this correlation breaks at high rotation rate with 
the possible appearance of a super active branch. Further 
at very high rotation rate saturation of the X-ray lumi- 
nosity seems to limit the validity of the scaling found at 
more moderate rotation rates (Pizzolato et al. 2003). For 
G type stars this saturation is found for rotation rate above 
35 kms -1 , for K type stars at about 10 kms -1 and for M 
dwarfs around 3-4 kms -1 (Browning 2008; Reiners, Basri 
& Browning 2009). How stellar magnetic flux scales with 
rotation rate is thus also important to understand since it 
is telling us how the magnetic field generated by dynamo 
action inside the stars emerges and imprints the stellar sur- 
face (Rempel 2008) and if it actually saturates. 

1.2. Possible theoretical interpretation 

Theoretical considerations to interpret stellar magnetism 
based on classical mean field a-uo dynamo models (Durney 
& Latour 1976, Knobloch et al. 1981; Baliunas et al. 1996) 
naturally find correlation between rotation rate and stel- 
lar activity. In particular it is found that both magnetic 
field generation and the dynamo number D (i.e. a Reynolds 
number characterizing the mean field a and uj dynamo ef- 
fects used in the models) vary with the rotation period 
of the star. This is due to the fact that in these models 
both effects are sensitive to the rotation rate of the star. 
The cj-effect is a direct measure of the differential rotation 



AQ established in the star. It is well known both theoret- 
ically and observationally that the differential rotation in 
the convective envelope of solar-type stars is directly con- 
nected to the star's rotation rate Qq (Donahue et al. 1996; 
Barnes et al. 2005; Ballot et al. 2007; Brown et al. 2008; 
Kiiker & Riidiger 2008). However, the exact scaling n r (i.e. 
AQ oc ^Q r ) is still a matter of debate among the observers 
and theoreticians, being sensitive to the observational tech- 
niques used and to the modelling approach. 

The a-effect used in mean field theory is related to he- 
lical turbulence and represents a measure or parameteriza- 
tion of the mean electromotive force (emf) (Moffatt 1978; 
Pouquet et al. 1976). It is thus also linked to the rotation 
rate of the star and the amount of kinetic helicity present 
in its convective envelope. However in Babcock-Leighton 
flux transport dynamo models the standard a-effect is re- 
placed by a surface term linked to the tilt of the active 
regions with respect to the east-west direction (i.e. Joy's 
law, Kosovichev & Stenflo 2008). This tilt is thought to be 
due to the action of the Coriolis force during the rise of the 
toroidal structures that emerge as active regions (D'Silva 
& Choudhuri 1993). Recent 3-D simulations in spherical 
shells (Fan 2008; Jouve & Brun 2007b, 2009) seem to in- 
dicate that this is not the only effect responsible for the 
observed tilt and that the twist and arching of the toroidal 
structures as well as the continuous action of the surface 
convection during the emergence have some influence on 
the resulting tilt. Thus the dependence with rotation rate 
of the Babcock-Leighton surface source term used in this 
class of mean field dynamo models is not straightforward 
to assess and needs to be studied in detail. 

Another important ingredient in flux transport models 
is the large scale meridional circulation (MC) used to con- 
nect the surface source term generating the poloidal field 
to the region of strong shear at the base of the convec- 
tion zone (i.e. the tachocline) where it will be subsequently 
sheared by the a;-effect in order to close the global dynamo 
loop (i.e. Bp i — > B tor — > Bpoi) . The meridional flow (or 
"conveyor belt") thus plays an important role in setting 
the cycle period of the global dynamo in this class of flux 
transport model. As a direct consequence it is natural to 
ask how the meridional circulation amplitude and profile 
change with the rotation rate and how these may influence 
the magnetic cycle period. 

Early work by Dikpati et al. (2001) have assumed that 
the amplitude of the meridional flow is proportional to the 
rotation rate of the star. With such scaling they seem to be 
able to reproduce some of the rotation vs cycle period cor- 
relation observed in solar type stars. Charbonneau & Saar 
(2001) have also studied the sensitivity of various mean field 
dynamo models in reproducing stellar activity observations 
assuming for the Babcock-Leighton type that the merid- 
ional flow amplitude increases either as Qq or as log(fio)- 
Similarly Nandy (2004) and Nandy & Martens (2007) have 
explored the solar-stellar connection with B-L mean field 
dynamo models and have reached the same conclusion: only 
a positive scaling of the amplitude of meridional flows with 
the rotation rate can reconcile the models with observations 
of magnetism of solar-like stars. However recent theoretical 
work by Ballot et al. (2007) and Brown et al. (2008) in- 
dicate that this is unlikely to be the right scaling. Instead 
3-D simulations reveal that the amplitude of the meridional 
flow weakens as the rotation rate is increased. 
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In this paper we wish to revisit the solar-stellar con- 
nection in the light of the latest theoretical and observa- 
tional work done in the community. We have thus designed 
a series of Babcock-Leighton models with the STELEM 
code (Emonet & Charbonneau 1998, private communica- 
tion; Jouve & Brun 2007a) that incorporate recent 3-D 
MHD results. By doing so we will study how these mean 
field models, quite successful at modelling the solar cycle, 
perform to explain the observed magnetic activity, cycle 
and field topology of solar type stars at various rotation 
rates. In §2 we discuss the 2-D mean field models used in 
this study and the way we incorporate 3-D results obtained 
with the ASH code (Clune et al. 1999; Brun et al. 2004). In 
§3 we present our series of stellar mean field dynamo mod- 
els, while in §4 we vary the input parameters to evaluate 
the robustness of our results. In §5 we discuss and put in 
perspective our results. 

2. The Models and Method 

2.1. The model equation 

To model the solar global dynamo operating in solar like 
stars, we use the hydromagnetic induction equation, gov- 
erning the evolution of the magnetic field B in response to 
advection by a flow field V and resistive dissipation: 



- = Vx(VxB) 



V x (r?V x B) 



(1) 



where 77 is the magnetic diffusivity. 

As we are working in the framework of mean-field the- 
ory, we express both magnetic and velocity fields as a 
sum of large-scale (that will correspond to mean field) 
and small-scale (associated with fluid turbulence) contri- 
butions. Averaging over some suitably chosen intermediate 
scale makes it possible to write two distinct induction equa- 
tions for the mean and the fluctuating parts of the magnetic 
field. The mean-field equation reads: 



9(B) 
dt 



V x ((V) x <B» + V x (v' x b') 

V x (77V x (B)) 



(2) 



where (B) and (V) refer to the mean parts of the mag- 
netic and velocity fields and v' and b' to the fluctuating 
components. A closure relation is then used to express the 
mean electromotive force (v' xb') in terms of mean mag- 
netic field, leading to a simplified mean-field equation. In 
this work we will replace the emf by a surface Babcock- 
Leighton term (Babcock 1961; Leighton 1969; Wang et al. 
1991; Dikpati & Charbonneau 1999; Jouve & Brun 2007a) 
as described in details below. 

Working in spherical coordinates and under the assump- 
tion of axisymmetry, we write the total mean magnetic field 
B and velocity field V as (for the sake of simplicity we omit 
( ) in the rest of the paper): 



B(r, M)=Vx (A (r, 0, f)e ) + S (r, 0, t)e 4 
V(r, 0) = v p (r, 0) + r sin fi(r, 0)e , 



(3) 
(4) 



Note that our velocity field is time-independent since 
we will not assume any fluctuations in time of the dif- 
ferential rotation Q or of the meridional circulation v n . 



Reintroducing this poloidal/toroidal decomposition of the 
field in the mean induction equation, we get two coupled 
partial differential equations, one involving the poloidal po- 
tential and the other concerning the toroidal field B^. 

cM <* ^(v 2 -i-)A - J R e ^-V(^)+C s 5(r,^^)(5) 



dt 



dt 



Vt 



Vt 



ZD 



(V 2 



1 d{zuB^) div/rjt) 



zu 



dr 



dr 



,B 



R e wv p -V(^)-R e B^V-v p 



(6) 



where w = rsin0, rjt is the turbulent magnetic diffu- 
sivity (diffusivity in the convective zone) and £(r, 0, B<j>) 
the Babcock-Leighton surface source term for poloidal field. 
In order to write these equations in a dimensionless form, 
we choose as length scale the solar radius (i?©) and as 
time scale the diffusion time (R%/rjt) based on the enve- 
lope diffusivity rjt = 10 11 cm 2 s _1 . This leads to the ap- 
pearance of three control (Reynolds numbers) parameters 
Cq = QoR^/rjt, C s = s Ro/rj t and R e = v R G /rjt where 
^o,5o, are respectively the rotation rate and the typical 
amplitude of the surface source term and of the meridional 
flow. 

Equations [5] and [6] are solved with the STELEM code 
(Emonet & Charbonneau 1998, private communication) in 
an annular meridional plane with the colatitude 6 G [0, tt] 
and the radius (in dimensionless units) r G [0.6, 1] i.e. from 
slightly below the tachocline (r = 0.7) up to the solar sur- 
face. The STELEM code has been thoroughly tested and 
validated thanks to an international mean field dynamo 
benchmark involving 8 different codes (Jouve et al. 2008). 
At 6 = and 6 = tt boundaries, both and B^ are set to 

0. Both A^ and B^ are set to at r = 0.6. At the upper 
boundary, we smoothly match our solution to an external 
potential field, i.e. we have vacuum for r > 1. As initial con- 
ditions we are setting a confined dipolar field configuration, 

1. e. the poloidal field is set to sin 0/r 2 in the convective zone 
and to below the tachocline whereas the toroidal field is 
set to everywhere. 



2.2. The physical ingredients 

The rotation profile used in the series of models dis- 
cussed in this work captures some realistic aspects of the 
Sun's angular velocity, deduced from helioseismic inversions 
(Thompson et al. 2003). We are thus assuming a solid ro- 
tation below 0.66 and a differential rotation above the in- 
terface (cf. top panel of Figure [T|). 



n(r,0) = ft c + I[l + er /(2- 



di 



(Q*Eq + &2 COS 2 + (24 COS 4 — Q c 



(7) 



with Q Eq = 1, tt c = 0.93944, r c 0.7, d x 0.05, 
a 2 = -0.136076 and a 4 = -0.145713. With this profile, the 
radial shear is maximal at the tachocline. As already dis- 
cussed in the introduction both observations (Donahue et 
al. 1995; Barnes et al. 2005) and 3-D numerical simulations 
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(Ballot et al. 2007; Brown et al. 2008) indicate that the in- 
ner rotation profile is likely to vary with the rotation rate. 
In the models discussed in section 3 we will keep the profile 
of differential rotation fixed, varying only the rotation rate 
Qo of the models. We will study the influence of varying 
the differential rotation profile in section 4.1. 



x [1 + ( ^% t] f}' 1 cos sin fll? (r c , fl, t) 




Fig. 1. Top panel: differential rotation profile used in the 
models GO. 25 to G10 discussed section 3. We zoomed on 
the tachocline region between r = 0.6 and r = 0.8. Bottom 
panel: meridional circulation used in the models GO. 25 to 
G10 discussed section 3. Plain (dashed) lines indicate coun- 
terclockwise (clockwise) circulation. 



We assume that the diffusivity in the envelope r] is dom- 
inated by its turbulent contribution whereas in the stable 
interior r] c <C r} t . We smoothly match the two different con- 
stant values with an error function which enables us to 
quickly and continuously transit from r] c to rjt i.e. 



rj(r) = rye + - (rft - rj c ) 



1 + erf 



(8) 



with f] c = 10 9 cm 2 s _1 and d = 0.03. 

In Babcock-Leighton (BL) flux-transport dynamo mod- 
els, the poloidal field owes its origin to the tilt of magnetic 
loops emerging at the solar surface. Thus, the source has to 
be confined to a thin layer just below the surface and since 
the process is fundamentally non-local, the source term de- 
pends on the variation of at the base of the convection 
zone (CZ). The expression is 



where r 2 = 0.95, d 2 = 0.01, B = 10 5 G. 

In these BL flux-transport dynamo models, meridional 
circulation is used to link the two sources of the magnetic 
field namely the base of the CZ and the solar surface. The 
structure of the meridional circulations in the Sun is less 
constrained than the internal differential rotation, as they 
possess slower flows that are difficult to detect helioseismi- 
cally. The nature of these flows in other stars is even less 
constrained. It is generally believed that in the Sun, the 
meridional circulations are likely to be composed largely of 
a single cell in each hemisphere, and recent solar simula- 
tions (Miesch et al. 2008) find such single-celled structures 
when the global-scale flows are averaged over long inter- 
vals but multi-cellular profiles can also be realized. In more 
rapidly rotating suns, 3-D models presently indicate that 
the circulations are likely to be multi-celled in both radius 
and latitude (Ballot et al. 2007; Brown et al. 2008). In the 
series of models discussed in section 3 of this paper, we re- 
strict ourselves to Babcock-Leighton flux-transport models 
that have a large single cell per hemisphere. We choose this 
simple profile because it is widely used in the community, 
although some studies have considered more complex pro- 
files (Bonanno et al. 2005; Jouve & Brun 2007). We are de- 
ferring to section 4.2 the analysis of complex multi-cellular 
profiles. As in the Sun, the meridional circulations are di- 
rected poleward at the surface and here they vanish at the 
bottom boundary (r = 0.6) and thus penetrate a little be- 
neath the tachocline into the radiative interior (cf. bottom 
panel of Figure [1]). To model the single cell meridional cir- 
culation we consider a stream function with the following 
expression: 



/ ( a\ 2 ( r ~ r b) 2 . ( r-r b 
ip(r, 0) — - sin ( 7T Z I cos V sin t 



7T (1 - r b ) 



l-n 



(9) 



which gives, through the relation v p = V x (^e^), the fol- 
lowing components of the meridional flow 



v r 



2(1 - rb) (r - r b ) 2 
7rr (1 — rb) 2 
r - r b 



X Sin 7T 



l-n, 



(3cos^6>- 1), 



(10) 



r 3r — rb . / r 

sin 7T— 

1 - r b V 1 

r7r (r — rb) 



rb 



1 - r b (1 - r b ) 
2(1 - r b ) (r - r b ) 
7rr (1 — rb) 



COS 7T 



1 - n> 



cos t) sin f 



(ii) 



with rb = 0.6. 
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3. Flux transport models as a function of rotation 
rate 

3.1. Previous results 

Jouve & Brim (2007a) have done a systematic study of the 
influence on the cycle period of the various parameters of 
Babcock-Leighton flux transport models. For the reference 
model possessing only one cell per hemisphere they con- 
firm the results of Dikpati & Charbonneau (1999) that the 
cycle period is strongly dependent on the meridional circu- 
lation amplitude. A least square fit through the ensemble of 
models they computed leads to the following dependency: 



-0.91, 



eye 



-0.013 - 

o V 



0.075 



-0.014 



(12) 



We can thus expect that as we increase the rotation rate 
the cycle period will shorten if all the other parameters are 
kept constant but unfortunately not as fast as what is ob- 
served (Noyes et al. 1984; Saar 2002). However we may also 
expect the meridional circulation to depend on the rotation 
rate, and the sign of the variation of with Qq will most 
certainly dominate over the weak dependence of P cyc on Qo. 
It is thus crucial to assess how the meridional circulation 
amplitude and profile vary with the rotation rate ^o- We 
will consider in the following sections both various ampli- 
tude for the meridional circulation at a given rotation rate 
and more complex profiles with several cells per hemisphere 
as was done in Jouve & Brun (2007a), to quantify how this 
may affect the stellar activity cycle. 

3.2. Implementing results from 3D simulations of rapidly 
rotating stars in our mean field dynamo models 

In order to get some physical insight on how the merid- 
ional circulation varies with the rotation rate we rely on 
recent 3-D numerical simulations of solar convection per- 
formed with the ASH code (Brun & Toomre 2002, Ballot 
et al. 2007, Brown et al. 2008). Indeed there are no reli- 
able observations of the meridional circulation in stars ex- 
cept for the Sun near the surface. The theoretical work of 
Brown et al. (2008), that studies the influence of the ro- 
tation rate on solar like stars (i.e. with the same thickness 
for the convection zone as our current Sun), will provide 
the required scaling. In this work several 3-D purely hydro- 
dynamical models have been computed at various rotation 
rates (from 1 to 10 times solar) and the redistribution of an- 
gular momentum, heat and energy assessed systematically 
to explain the differential rotation and meridional circu- 
lations achieved in the simulations. Figures 10 and 12 of 
Brown et al. (2008) are particularly relevant to our study 
and we have used the data from all the models reported 
in these figures plus more slowly rotating cases (i.e. with a 
rotation rate lower than the solar rate) that are currently 
being studied (Brown 2009). 

We have thus scaled the meridional circulation by taking 
into account the way meridional circulation kinetic energy 
(MCKE) is found to change with rotation in the published 
3-D models of Brown et al. (2008). We are confident that 
this scaling for MCKE still holds in the magnetized version 
since recent 3-D models of dynamo action in rapidly ro- 
tating suns find that the weak meridional circulations are 
essentially unmodified by the strong magnetic fields that 
are generated (Brown et al. 2009). In contrast, the differen- 
tial rotation appears to be substantially different from that 





On /Oo 


MCKE/10 4 


scale 


"R /If) 2 


n n /i n 5 


GO. 25 


0.25 


23.4 


28.7 


20.1 


0.351 


GO. 5 


0.5 


16.2 


23.9 


16.7 


0.702 


GO. 6 


0.6 


12.3 


20.9 


14.6 


0.842 


GO. 75 


0.75 


4.75 


12.9 


9.06 


1.05 


GO. 9 


0.9 


3.09 


10.4 


7.29 


1.26 


Gl 


1 


2.84 


10.0 


7.00 


1.40 


G1.25 


1.25 


2.18 


8.76 


6.13 


1.75 


G1.5 


1.5 


1.98 


8.36 


5.85 


2.11 


G1.75 


1.75 


1.68 


7.70 


5.39 


2.46 


n 9 

V I Z, 


2 


1.52 


7 ^ 

I .00 


O. J-O 


2.81 


G3 


3 


1.11 


6.25 


4.38 


4.21 


G4 


4 


0.767 


5.20 


3.64 


5.62 


G5 


5 


0.711 


5.01 


3.51 


7.02 


G6 


6 


0.626 


4.70 


3.29 


8.42 


G7 


7 


0.497 


4.18 


2.93 


9.83 


G10 


10 


0.359 


3.56 


2.49 


14.0 



Table 1. Summary of the different cases studied and the 
values of R e and Cq used in the 2D Babcock-Leighton mod- 
els. The value of rjt is fixed to 10 11 cm 2 .s _1 and the value of 
so is fixed to 50 cm.s -1 . The stellar rotation rate is given 
in terms of the solar rotation rate Q$- v s e cale represents the 
amplitude of vq at the surface at 45° and is imposed equal 
to 10 m.s -1 for case Gl, representing the reference solar 
mean field dynamo model. 



found in hydrodynamic cases, and its scaling with rotation 
rate is sensitive to magnetic fields that are realized. We will 
thus not use the scaling of with found in Brown et 
al. (2008) when we will vary the differential rotation profile 
in our models but instead rely on observations (see section 
4.1 below). 

In Table [U we list the solar-like mean field dynamo mod- 
els that we have computed at various rotation rate along 
with the associated Reynolds numbers R e and Cq that we 
have used in our 2-D models. A quick look at the table 
already reveals that the meridional circulation amplitude 
decreases with rotation rate. This may come as a surprise 
as one could have expected faster meridional flows as Qo 
is increased. However a careful study, via the longitudinal 
component of the vorticity equation, of the maintenance 
of the meridional circulation (which arises from small un- 
balance between large forces; see for instance Miesch 2005; 
Brun & Rempel 2008), reveals that it actually weakens with 
faster rotation rate. It is thus likely given the strong neg- 
ative dependency of P cyc with that as we increase the 
rotation rate, the cycle period will increase rather than de- 
crease as observed. However it is not obvious to anticipate 
the resulting field topology of the various models without 
actually computing them. 

Incorporating the values of the dynamo numbers R e and 
Cq for the different stars in our flux-transport BL dynamo 
models enables us to study the behaviour of the mean mag- 
netic field in such stars in the framework of mean-field the- 
ory. All the dynamo models we have studied give cyclic 
behaviours with different periods and different intensity of 
the mean toroidal and poloidal magnetic fields. We wish to 
stress that even though the naming scheme of the models 
discussed in this paper is similar to that used in Brown et 
al. (2008), the models are distinct since the ones discussed 
in this work are 2-D not 3-D. 
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Fig. 2. Butterfly diagrams for cases GO. 5, Gl, G5 and G10. 
The top panel shows the time-latitude cut of the toroidal 
field at the base of the convection zone and the bottom one 
represents the evolution of the radial field at the surface. 
Time is in years. The scale of the color table is the same 
for all: between —5 x 10 3 G and 5 x 10 3 G for the radial 
field and between —5 x 10 5 G and 5 x 10 5 G for the toroidal 
field. 



Figure [2] gives examples of 4 different butterfly diagrams 
obtained for more or less rapidly rotating stars. This figure 





Fig. 3. Representative solution of a BL model applied to 
a rapidly rotating star. This is case G5, with a rotation 
rate of 5 Vts • We show the evolution of the toroidal 
field (filled colour contours, with polarity indicated: red, 
positive; blue, negative) and of the poloidal field (contour 
lines) over a period of 20.8 years, or half a magnetic cycle 
in this simulation. A potential field solution to the poloidal 
field is extrapolated from the surface to 1.2R Q . 



shows the evolution in time of the toroidal magnetic field 
at the bottom of the CZ and of the radial field at the top of 
the domain. We note that the behaviour of the field is very 
similar in all these models, except that the period is clearly 
varying as well as the intensity of the magnetic field (the 
scale of the color table is the same for all panels). Indeed, we 
note that the magnetic cycle period is increased when the 
rotation rate is increased. This is a direct consequence of the 
fact that the cycle period is almost inversely proportional 
to the amplitude of the meridional flow in BL models. 

In order to illustrate the temporal evolution of the mag- 
netic field in the meridional plane, we display on Figure 
[3] a temporal sequence over one starspot cycle for model 
G5. The advection of the field by the single-celled merid- 
ional circulation is evident confirming that our dynamo 
model is dominated by advection and not by diffusion. We 
wish to stress here that our particular choice of advection- 
dominated regime is not the only possibility at hand to 
model the transport occurring in the solar convection zone 
and the manner that the models link the Babcock-Leighton 
surface source term to the deep tachocline. For instance 
Chatterjee et al. (2004), Kapyla et al. (2006), Yeates et al. 
(2008), Guerrero & Gouveia Dal Pino (2008) have studied 
how transporting the poloidal field from the surface to the 
tachocline via turbulent diffusion or turbulent pumping in- 
fluence the butterfly diagram, the polarity of the field and 
the cycle period. We intend to look at these solutions in the 
near future. 

Returning to Figure [3] we also note that the poloidal 
field tends to be more concentrated near the polar region 
in agreement with the stronger polar branch seen in the cor- 
responding butterfly diagram. Indeed in Figure [2] we clearly 
note that the magnetic intensity is increased when the ro- 
tation rate is increased, which is in agreement with the 
observations. We finally also note on the poloidal field that 
the contrast between the polar intensity and the intensity 
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close to the equator is diminished when the rotation rate is 
increased. 

3.3. Field topology and poloidal vs toroidal field ratio 

It is clear that the evolution of the field topology as a func- 
tion of the rotation rate is of great importance. We are thus 
summarizing in Table [2] and on the curves of figure [H how 
the field intensity and topology is changing with and 
the meridional circulation amplitude. We show for instance 
the variation of the magnetic cycle period with respect to 
the amplitude of the MC, with respect to the rotation rate 
and the variation of the maximum of toroidal field at the 
base of the CZ with respect to the rotation rate. 



Case 


n /n s 


Cycle period 


Max Bcf, 


Bpol/Btor 


G0.25 


0.25 


4.65 


1.56 


0.398 


G0.5 


0.5 


5.42 


3.66 


0.179 


G0.75 


0.75 


9.37 


6.00 


0.0889 


Gl 


1 


11.7 


7.50 


0.0555 


G1.5 


1.5 


13.6 


9.50 


0.0319 


G2 


2 


15.1 


11.0 


0.0222 


G3 


3 


17.3 


14.3 


0.0120 


G5 


5 


20.8 


20.5 


0.00500 


G7 


7 


24.5 


26.6 


0.00380 


G10 


10 


28.5 


34.4 


0.00210 



Table 2. Starspot cycle period (in year), maximum of 
toroidal field (in 10 4 G) at the base of the CZ and ratio 
between poloidal and toroidal field for 10 different cases 
with various rotation rates. 



More precisely on the first (top) panel of Figure [H we 
recover the almost proportional relationship between the 
magnetic cycle period and the inverse of the meridional flow 
amplitude. This again shows the major impact of the merid- 
ional circulation in these Babcock-Leighton dynamo models 
on the cycle period. As stated in Brown et al (2008), the 
meridional circulation kinetic energy in 3-D hydrodynami- 
cal (HD) models of rapidly rotating stars seems to scale as 
approximately ^q 0,9 , leading the averaged meridional flow 
amplitude to scale as Qq 0A5 . Indeed, on the second panel 
of Figure [4] where we show the cycle period plotted against 
the rotation rate of the star, we naturally recover a scaling 
in ^q- 45 . A star rotating at SCls for example would have a 
cycle period of about 26 years approximately 4 45 = 1.87 
times the period of 15.1 years of case G2 rotating at twice 
the solar rate. 

More interestingly, if we now focus on the third and 
fourth (bottom) panels of Figure 0J we note that the ampli- 
tude of the toroidal field possesses an almost linear depen- 
dency on the rotation rate ^o- This property is due mainly 
to two effects: a) a larger o;-effect (via its dependence on 
the Reynolds number Cq) in more rapidly rotating stars 
and b) the fact that the reduction of the meridional flow 
amplitude allows the poloidal field to be stretched and am- 
plified over a longer period of time (see Nandy (2004) and 
Yeates et al. (2008) who are finding similar effects). This 
is evident by comparing on the third panel of Figure [H the 
solid and dashed curves. The difference between the two 
curves is the amplitude of the meridional flow as a function 
of rotation rate that we kept constant for the cases shown 
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Fig. 4. Top panel: Variation of the magnetic cycle period 
wih respect to the amplitude of the MC and (2nd panel) to 
the rotation rate for the cases Gl to G10. On the third panel 
we display the variation of the maximum of toroidal field at 
the base of the CZ for two series of models. The solid line 
corresponds to the models Gl to G10 using the scaled down 
meridional flow as a function of the rotation rate whereas 
the dotted line corresponds to another set of models for 
which was kept constant. On the bottom panel, we show 
the ratio Bpol/Btor wih respect to the rotation rate for 
cases Gl to G10. 
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using dashed line (these models were used to determine the 
least square fit of equation (12)). We clearly see that both 
series of models lead to larger toroidal field as we increase 
the rotation rate but for the series with the scaled down 
meridional flow, the strengthening of is even higher. As 
a result, in our set of 2-D BL dynamo models the poloidal to 
toroidal field ratio is found to increase almost linearly with 
the rotation period (inversely proportional to the rotation 
rate as seen on the last panel of figure [4|) . This is in rea- 
sonable agreement with the observations of rapidly rotating 
solar-like stars by Petit et al. (2008). However the scaling of 
Bpol/Btor found in recent 3-D MHD simulations computed 
by Brown et al. (2009) is not linear. Moreover, the obser- 
vations indicate that the energy stored in the large-scale 
toroidal component of the magnetic field tends to domi- 
nate for rotation periods lower than 12 days (correspond- 
ing approximately to a rotation rate of twice the solar one). 
In our cases, the toroidal field is always dominant, even for 
the case with 0.25^5 where the poloidal field hardly reaches 
40% of the toroidal field. 

This discrepancy could be due to the fact that the values 
quoted in Petit et al. (2008) are related to observationally 
detectable fields, i.e. to surface fields; whereas the values we 
quote in this work are related to the magnetic fields in the 
whole convection zone. If the strong toroidal fields stored 
at the base of the convection zone are not able to reach the 
solar surface (as it is the case under certain conditions, see 
Jouve & Brun 2009 for more details), the surface measure- 
ments could underestimate the amount of magnetic energy 
stored in the toroidal fields in the whole stellar convective 
envelope. 
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4. Varying the model parameters 

From the above set of models it is clear that just increas- 
ing the rotation rate of the model and/or implementing a 
scaled meridional circulation amplitude is not sufficient to 
reproduce the fast (close to linear) reduction of the stellar 
cycle period with the rotation rate. We thus wish in this 
section to modify all the other parameters in order to find 
the best set of models to reproduce the observations. 

The first easiest modification is to change the value of 
the diffusivity or of the constant sq used in the model. From 
the least square fit shown in equation ([T2]) the weak depen- 
dency of P cyc on r] or so indicates that it is unlikely that 
such a modification will be able to reproduce the observa- 
tions. We thus need to look for more sophisticated solutions. 
In order to make the discussion easy to follow and to re- 
duce the number of models considered we will now limit 
our study to cases rotating at five times the solar rate. Of 
course our results are easily transposable to other rotation 
rates. 

4.1. Modifying the differential rotation profile 

As discussed in the introduction it is clear that in low mass 
stars the differential rotation is expected to change with 
the rotation rate. Given the key role played by the o;-effect 
in mean field dynamo, we propose to compute a new model 
at 5 times the solar rate with a modified differential rota- 
tion, namely case G5om. But what scaling of AQ with Qq 
should we use? In their 3-D hydrodynamic models, Brown 
et al. (2008) find a scaling of A^ oc ^q" 3 . However due to 
the nonlinear feed back of the Lorentz force on the mean 



Fig. 5. Top panel: New differential rotation used in the 
model G5om zoomed on the portion between r = 0.6 and 
r — 0.8. For that case AQ has been increased by about a 
factor of 3 following Donahue et al. (1996) scaling. Bottom 
panel: Butterfly diagrams for that same case G5om with 
increased AQ. The format is the same as Figure [2j 



longitudinal flow, this scaling is not recovered in the MHD 
models computed in Brown et al. (2009). We thus choose 
here not to use the theoretical scaling for AVt of Brown 
et al. (2008) since it is somewhat uncertain and depends 
on whether or not dynamo action is realized in the simula- 
tions. The observational situation is not clearer with either 
a weak scaling close to zero (i.e. the differential rotation 
does not vary with the rotation rate; Barnes et al. 2005) 
or strong scaling as in Donahue et al (1996). Given the 
current debate going on among the observers on the ex- 
act variation of A^ with Qq we have some freedom in the 
choice of our variation. We choose here to follow the scaling 
proposed by Donahue et al. (1996), i.e. A^ oc Qq- 7 since it 
is the most extreme variation found among the observers 
and as a consequence it will maximize the effect of varying 
the differential rotation profile. The new differential rota- 
tion profile is shown on the top panel of Figure [5] and is 
directly comparable with the one shown on the top panel 
of Figure [TJ We also display on the bottom panel of Figure 
[5] the butterfly diagram for this new model G5om rotating 
at 5 times the solar rate with increased AQ. 

We see that the cycle period is almost not affected by 
the increase of the latitudinal shear when comparing the 
butterfly diagram of Figure [5] with the third panel of Figure 
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[21 This confirms that the o;-effect plays little role in setting 
the cycle period in BL dynamo models. On the contrary 
the ratio between Bpol and Btor is changed. Indeed, for this 
case, we go from a Bpol to Btor ratio of 0.05 in the standard 
case G5 to a ratio of 0.019 in G5om. We thus proceeded in 
this case to a significant increase of the toroidal component 
of the magnetic field, keeping the value for the poloidal field 
approximately the same. This is due to the fact that this 
new class of models is generating stronger Btor as can be 
expected with a larger shear. A precise observation of this 
ratio may thus be of utility to put constraints on stellar 
internal differential rotation. 

With a larger latitudinal shear it is expected that the 
tachocline at the base of the stellar convection zone will 
also be modified. The tachocline is the transition region be- 
tween the differentially rotating convective envelope of low 
mass star and their supposedly uniformly rotating radia- 
tive interior. How the tachocline thickness is varying with 
the rotation rate and the differential rotation of the star is 
not easy to determine. Based on the tachocline model of 
Spiegel and Zahn (1992) and the observations of Donahue 
et al. (1996) for the dependency of AQ with Qq, Brun et 
al. (1999) have shown that the tachocline thickness for a 
younger more rapidly rotating Sun could be thicker. They 

find a scaling for d oc ^q 1T±01 ^ 4 when taking into account 
all the dependencies. We have thus computed a new model 
G5omt, with both a larger tachocline (about twice as large 
in this case than in case G5om, following the scaling law 
quoted above) and stronger differential rotation (identical 
to case G5om shown on Fig. [5j with a AQ about three 
times the solar one). The result in term of the cycle period 
is again difficult to distinguish from the standard G5 model, 
that had only Cn and the meridional circulation amplitude 
changed. This confirms that the tachocline plays almost 
no role in Babcock-Leighton dynamo models in setting the 
cycle period. On the contrary the ratio Bpol/Btor is mod- 
ified being intermediate between the cases G5 and G5om. 
The value of the ratio is now of about 0.02, very close to 
the case where the differential rotation was increased, but 
with a slight reduction of the toroidal component, due to 
the softer transition (i.e. thicker tachocline). In summary 
neither models G5om nor G5omt possess a cycle period in 
agreement with the observations and we thus need to mod- 
ify another physical ingredient. 

In the standard Babcock-Leighton model we are thus 
left with two possibilities if we want to reconcile the models 
and the observations: modifying the profile of meridional 
circulation or introducing the effects of a faster rotation rate 
on the surface source term. Indeed, for example it has been 
shown by D'Silva & Choudhuri (1993) that the Coriolis 
force could strongly act on the tilt angle of emerging active 
regions, thus modifying the source term for poloidal field. 
This feature will have to be checked in future works, with 
the help of 3-D MHD models such as the ones of Jouve 
& Brun (2009). For this work, since the dependence of the 
cycle period on the amplitude of the meridional flow is very 
strong, we decide to focus on the first solution: modifying 
the MC profile. 

4.2. Varying the number of meridional circulation cells 

The strong relation between the meridional circulation 
speed and the cycle activity is an indication that shorten- 



ing the advection path may result in shorter cycle period. 
Dikpati et al. (2004), Bonanno et al. (2005) and Jouve & 
Brun (2007a) have studied in detail the influence of various 
meridional circulation profiles on the butterfly diagram and 
other large scale magnetic properties of the Sun. 

Stream function Latitudinal velocity 




latitude 

Fig. 6. Multicellular stream function for the case rotation 
at 5 times the solar rate G5mc20. Two cells with latitude 
have been used, with the first counterclockwise cell near 
the equator closing at 20 deg of latitude. On the left panel 
we show the streamfunction and on the right panel the 
latitudinal velocity at 3 depths. 

In particular Jouve & Brun (2007a) have shown that 
adding cells in radius slows down the cycle whereas as the 
two other groups also showed, adding cells in latitude ac- 
celerates it. If one seeks to shorten the cycle length as it is 
the case in this study, only adding cells in latitude will have 
the expected impact. We have thus computed a new set of 
simulations, rotating at five times the solar rate in which 
we have only modified the meridional circulation in latitude 
(leaving all the other parameters the same as in case G5) 
by considering two cells in latitude per hemisphere with 
various extension instead of just one cell. 

Following the 3-D results of Brown et al. (2008) and 
particularly figure 11 of their article, we choose to intro- 
duce a meridional circulation with a counterclockwise cell 
at lower latitudes and a clockwise cell at higher latitudes in 
the Northern hemisphere, this profile being antisymmetric 
with respect to the equator. Since uncertainties are signif- 
icant concerning the extension of the low-latitude merid- 
ional cell, we have computed several models at 5 times the 
solar rate including a low-latitude cell extending to 45°, 
40°, 30° and 20°, named respectively G5mc45, G5mc40, 
G5mc30 and G5mc20. For the sake of clarity we have de- 
cided to illustrate and discuss in length in this section only 
one case, namely case G5mc20, the model with the less 
extended lower meridional circulation cell. The meridional 
circulation profile for that case is shown on Figure [6] for 
illustration. In all models we have set the latitudinal veloc- 
ity such that it has the same maximal amplitude in the two 
counter-cells which correspond to that in the case with one 
single cell (i.e. case G5) keeping all the other parameters 
the same. 

When comparing the various multi-cellular stellar dy- 
namo simulations, we find that the scaling of the cycle pe- 
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Fig. 7. Butterfly diagram of the case G5mc20, rotating at 5 
times the solar rotation rate and possessing a multicellular 
meridional flow as shown on Figure [6j 

riod with the extension of the low-latitude cell is almost 
linear, the starspot cycle period being respectively of 9.9, 
9.2, 7.4 and 5.2 years for those different models. This is 
mainly due to the fact that the advection path is reduced 
by the presence of counter-cells and will shorten the dis- 
tance each component of the magnetic field has to travel 
before reaching the region where it will be turned into the 
other component. According to the previous studies and 
scaling laws of Noyes et al. (1984), Saar & Brandenburg 
(1999), Charbonneau & Saar (2001) and Saar (2002), the 
cycle period for a star rotating at 5 times the rotation rate 
should be around 3 years. We are getting very close to this 
value with a low-latitude cell extending to 20° in latitude 
as in case G5mc20. We now focus on butterfly diagram and 
the evolution of the field in the meridional plane of this 
model, as shown on Figures [71 and 

On the butterfly diagram shown on Figure [3 we clearly 
see that we have significantly reduced the length of the ac- 
tivity (starspot) cycle: the period is now 5.2 yr instead of 
22 yr. The cyclic activity tends now to be more concen- 
trated at low latitudes. The temporal evolution of the field 
seems to be less steady if we consider the variations of in- 
tensity which appear from one cycle to the next, visible on 
the butterfly diagram through changes in color between the 
different activity cycles. The radial field appears to be even 
more busy, exhibiting a periodic but yet complex pattern at 
low latitudes. We moreover note that above 55°, the cyclic 
activity is lost and we are left with some regions where the 
magnetic field evolves much more slowly than the time scale 
of the main cycle. 

Looking at Figure [8] will help us to have more insight on 
how the field is organized in the convective envelope and 
how it is advected by the new meridional circulation. Again, 
we clearly see that the evolution of the field is much more 
complicated than in the standard case. We can still follow 
the evolution of each component of the field on the various 
panels of the figure. Indeed we can distinguish the effect 
of the low-latitude cell which advects the newly generated 
poloidal field (counterclockwise on the first panel) down- 
wards while the toroidal field of the previous cycle (posi- 
tive i.e. represented by reddish colours) is being advected 
at the base of the convection zone to reach the equator and 




Fig. 8. Same as Figure [3] but for the case G5mc20 rotating 
at 5 times the rotation rate and possessing a multicellular 
meridional circulation as shown on Figure [6l The evolution 
is shown over a temporal interval of 5.2 years corresponding 
to half a magnetic cycle or one starspot cycle. 



replace the toroidal field of the previous cycle at the end 
of the starspot cycle (last panel). We note, in agreement 
with the butterfly diagram, that the advection at higher 
latitudes is much slower and as a consequence we do not 
clearly see any imprint of the activity cycle in the polar 
regions. This is due to the fact that the same amplitude 
of the flow is introduced for both cells and that the dis- 
tance traveled by the field at low latitudes is much shorter, 
thus leading to a difference in the time-scale for the advec- 
tion of the magnetic field in the two cells. We also wish 
to reemphasize that we did not introduce multi cells in ra- 
dius since Jouve & Brun (2007a) showed that it results in a 
slowing down of the cycle period rather than a speeding up 
as required in this study for fastly rotating stars. In sum- 
mary adding cell in latitude clearly helps shortening the 
cycle for fastly rotating stars. For the stars rotating slower 
than the Sun it may be necessary to add cells in radius if 
one wants to increase the cycle. Indeed 3-D models indicate 
that slower the star rotates faster its meridional circulation 
is and thus shorter its cyclic period will be as demonstrated 
by the butterfly diagram shown on the top panel of Figure 
IH for case G0.5. 

5. Summary and Conclusions 

The idea of this work was to better understand stellar mag- 
netism, using both 2-D and 3-D calculations. To do so, we 
have introduced into Babcock-Leighton (B-L) flux trans- 
port dynamo models, results from 3-D global hydrodynam- 
ical simulations by Brown et al. (2008). To be more specific 
we have used these simulations to determine the variation 
of the amplitude of the meridional circulation with rota- 
tion rate and complemented our study with observations 
of the variation of the angular velocity contrast AQ with 
rotation rate. Using these scaling laws we have built a se- 
ries of 2-D mean field stellar dynamo models in which we 
have modified in turn the rotation rate, the amplitude of 
the meridional circulation and of the angular velocity gra- 
dient. This has allowed us to test the validity of this type of 
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models on solar-like stars with faster rotations and intense 
magnetism. 

Our study shows that an increase in the rotation rate 
leads to a larger a;-effect which in turn modifies the ratio 
between poloidal and toroidal fields with a trend similar 
to that found observationally in a small sample of rapidly 
rotating suns (Petit et al. 2008). Other authors have also 
studied such trends using solar mean field B-L flux trans- 
port dynamo models (see for instance Charbonneau & Saar 
2001; Dikpati et al. 2001; Nandy & Martens 2007). They 
also find that the toroidal field grows in amplitude with the 
star's rotation rate. The difference between our work and 
previous studies comes from the choice made in varying the 
meridional circulation amplitude with rotation rate. In 3-D 
global calculations an increase of the rotation rate leads 
to a decrease in the meridional flow amplitude whereas 
in previous studies the opposite was assumed since such 
knowledge was not yet available. Given that the cycle pe- 
riod in advection-dominated B-L dynamo models is almost 
inversely proportional to the intensity of the meridional cir- 
culation, we find that our models predict slower magnetic 
cycles for more rapidly rotating stars. This is not consistent 
with the observational data collected at Mount Wilson since 
the 1960's (e.g. Noyes et al. 1984, Baliunas et al. 1995). It 
thus seems difficult to reconcile advection-dominated B-L 
models widely used in the solar dynamo community with 
stellar activity observations if we use scaling laws relating 
rotation rates and meridional flow amplitudes given by con- 
sistent 3-D models. Of course the systematic observation of 
the meridional circulation realized in solar-like stars rotat- 
ing at various rotation rate would be of great use and could 
solve the apparent discrepancy if we indeed found that such 
flows are not so sensitive to the rotation rate. Unfortunately 
such observations are not yet available. Hopefully with the 
Corot and Kepler satellites now in orbit this situation may 
change in the near future. 

We have thus tried to investigate ways to improve the 
agreement between 2-D models and observations in refining 
our calculations. A least-square fit on the various param- 
eters of the standard model shows that we cannot rely on 
a modifications of fio, so and r] to act on the cycle period, 
as the latter is almost completely determined by the am- 
plitude of the meridional flow (see Equation 12). We thus 
computed more sophisticated models where the differential 
rotation and meridional circulation profiles were modified. 

First, we computed a model at 5 times the solar rota- 
tion rate where the modified contrast in the angular velocity 
was taken into account. Some uncertainties both from an 
observational or a theoretical standpoint still exist concern- 
ing the scaling between AQ and Qq but we find that the 
effect on the cycle period is negligible, even when we con- 
sider the most extreme observation of the variation of the 
differential rotation contrast with rotation rate (Donahue 
et al. (1996)). Modifying the tachocline thickness as a func- 
tion of rotation rate according to the scaling law of Brun et 
al. (1999) does not help to improve the relationship between 
the rotation rate of the star and the cycle period either. We 
show that the main impact of a modification of the profile 
of differential rotation is to decrease the percentage of the 
poloidal component in the total magnetic field again be- 
cause of the increased efficiency of the source for toroidal 
field: the o;-effect. The modification of the source term for 
poloidal field (decay of active regions in Babcock-Leighton 
models and the kinetic helicity of the flow in a-uj models) 



due to an increase of the rotation rate still needs to be ad- 
dressed to study the evolution of the poloidal to toroidal 
field ratio in those refined cases. 

Given the strong dependency of the magnetic cycle pe- 
riod on the meridional circulation and the fact that 3-D 
models exhibit multicellular patterns of the flow in lati- 
tude, it is reasonable to think that modifying the profile 
of the MC according to 3-D models may help to get closer 
to the observational data. Indeed, we show that adding a 
counter cell in latitude which extends until about 20° in 
both hemispheres reduces drastically the cycle period. We 
find a quasi-linear relationship between the extension of the 
counter cell and the cycle period, which is very promising 
to reconcile Babcock-Leighton models with observations of 
stellar activity. Another solution could be to change the 
way poloidal magnetic fields are transported from the sur- 
face down to the tachocline. In Yeates et al. (2008) the 
transport is dominated by diffusion to link the distant mag- 
netic field generation regions. In their case the cycle period 
seems less sensitive to the meridional flow amplitude and 
profile assumed and depends more on the diffusion profile. 
Guerrero & Gouveia Dal Pino (2008) have introduced in 
Babcock-Leighton dynamo equations an extra term repre- 
senting the pumping of magnetic fields due to turbulent 
convection (inspired by the work of Kapyla et al. 2006). 
This new process leads to vertical and horizontal transport 
of the fields. Their model seems to show a significant de- 
crease of the cycle period dependency on the meridional 
flow amplitude (from Vq ' 9 to Vq ' 12 , although only a shal- 
low meridional flow is considered, which may influence the 
scaling law) and that the turbulent pumping is now respon- 
sible for regulating the cyclic activity. Reintroducing the 
meridional flow amplitude scalings deduced from 3-D cal- 
culations in this type of models would thus have much less 
impact on the cycle period. The changes in the pumping 
coefficients with respect to the rotation rate would have to 
be investigated to see if some of the observations of stellar 
magnetic activity could be recovered. 

In conclusion, advection-dominated Babcock-Leighton 
flux-transport models may offer insights into the dynamo 
action occurring within the convection zones of solar-type 
stars. However, it is necessary to make serious modifications 
to these models if we wish to reconcile their predictions with 
observations of real stellar magnetism. 
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